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1.  Introduction 


Brillouin  scattering  is  one  of  the  lowest  order  nonlinear  effects  that  arises  in  optical  fibers  and  it, 
thus,  limits  the  transmitted  laser  power  PB.  Above  a  threshold  laser  power,  Pt h,  the  Brillouin 
scattering  becomes  stimulated  rather  than  spontaneous,  and  the  Stokes  backscattered  power  Ps 
rises  dramatically.  Various  methods  have  been  devised  to  increase  the  threshold.  In  the  present 
report,  we  analyze  the  effect  of  linearly  sweeping  the  center  frequency  of  the  pump  laser  (while 
keeping  its  linewidth  narrow)  at  a  fast  chirp  rate  (3  (in  Hz/s).  To  be  effective,  the  chirp  must  be 
large  enough  that  the  pump  laser  gets  swept  out  of  the  Brillouin  gain  bandwidth  Avq  within  the 
transit  time  of  the  fiber,  tf  =  nL/c,  where  L  is  the  length  of  the  fiber,  n  is  the  core  refractive  index, 

and  c  is  the  speed  of  light.  That  is,  one  requires  (i  to  be  large  enough  that  q  »  q,  where  tc  = 
AvB/2 /?.  Here  we  analyze  values  of  ft  up  to  1016  Hz/s,  corresponding  to  the  upper  range  of 

chirping  that  has  been  demonstrated  to  date  for  a  diode  laser  by  ramping  its  current. 

Section  2  explains  how  Brillouin  scattering  is  simulated  numerically  in  a  fiber  pumped  by  a 
chirped  laser.  Section  3  then  presents  the  simulation  results  for  a  6-km  fiber.  Such  long  fibers 
have  a  low  threshold  because  Pth  ~  2  lA/goL  in  the  unchirped  case,  where  A  is  the  modal  area  and 
g0  is  the  peak  Brillouin  gain  coefficient.  Consequently,  low  laser  powers  can  be  used  to 

experimentally  verify  the  simulations.  Good  agreement  is  found  between  theory  and  experiment. 
Next,  section  4  presents  results  of  the  time -dependent  simulations  for  a  short  fiber  ( L  =  17.5  m) 
having  characteristics  similar  to  those  used  for  high-power  laser  delivery.  Correspondingly 
larger  chirps  are  needed  to  suppress  stimulated  Brillouin  scattering  (SBS)  in  this  case.  The 
results  are  in  good  agreement  with  a  simpler  adiabatic  model  for  the  range  of  chirps  explored. 


2.  Theory  for  the  Dynamic  Simulations 


The  Brillouin  scattering  is  modeled  by  three  complex  coupled  partial  differential  equations 
(PDEs): 
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Here,  the  laser  electric  field  EL  ( z,t ),  Brillouin  Stokes-shifted  electric  field  Es  ( z,t ),  and  density 
variation  p(z,t)  of  the  fiber  from  its  mean  value  po  depend  on  time  1  and  longitudinal  position  z 
(varying  from  z  =  0  at  the  front  face  of  the  fiber  to  z  =  E  at  the  rear  face).  The  spectral  full- 
width-at-half-maximum  (FWHM)  of  the  spontaneous  Brillouin  peak  is  AvB.  For  a  silica  fiber  at 
an  incident  laser  wavelength  of  Al  =  1.55  pm,  the  index  is  n  =  1.447,  the  mean  density  is  /p  = 
2210  kg/m  ,  and  the  speed  of  sound  is  o  =  5960  m/s  (i).  The  loss  coefficient  in  the  fiber  is  taken 
to  be  a  =  0.2  dB/km  =  0.0461/km.  The  optic  coupling  parameter  is 
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where  the  electrostriction  coefficient  for  a  silica  optical  fiber  is  y=  0.902  (2)  and  the  polarization 
is  presumed  to  be  completely  scrambled  in  the  fiber  so  that  M  =  1.5  (2).  The  acoustic  coupling 
parameter  is 
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where  Sq  is  the  permittivity  of  free  space.  The  Langevin  noise  source/ (z,/)  is  delta-correlated  in 
time  and  space  such  that  ( 4 ) 

(/  (z,  t)f  *  (z\  t'))=Q  8{z  -  z'W  -  f)  (4) 

where  the  thermal  phonons  are  described  by  the  strength  parameter 

AMp oAvb  (5) 
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Here,  k  is  Boltzmann’s  constant,  T  =  293  K  is  room  temperature,  and  A  is  the  fiber  modal  area. 


The  three  PDEs  in  equation  1  are  solved  by  iterated  finite-difference  approximations  on  a  grid  of 
time  and  space  points  (5).  The  spatial  step  size  dz  is  chosen  to  be  small  enough  that  a  factor-of- 
two  reduction  does  not  change  the  final  values  in  equation  7  by  more  than  10%.  The  temporal 
step  size  is  dt  =  n  dzJc.  The  boundary  conditions  are  Es(L,t )  =  0  and 
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where  Po  is  the  constant  incident  laser  power  at  z  =  0.  The  imaginary  part  of  the  argument  of  the 
exponential  gives  the  chirped  phase,  equal  to  the  time  integral  of  the  change  in  the  laser 
frequency  at  a  rate  of  / 3  in  Hz/s,  relative  to  the  frequency  at  the  rear  face  of  the  fiber  at  t  =  0.  The 
three  complex  fields  are  first  calculated  across  the  fiber  up  to  t  =  20/ 1  to  ensure  that  the  initial 
relaxation  oscillations  (6)  have  decayed  away.  The  equations  are  then  iterated  over  five  more 
transit  times  to  acquire  a  statistical  average.  The  transmitted  laser  power  PL  and  reflected  Stokes 
power  Ps  are  computed  by  averaging  over  those  additional  steps, 
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PL  =  \nc£QA{\EL{L,t)\2)  and  Ps  =  ±nc£0A(\Es(0,t)\2).  (7) 

The  results  were  checked  by  verifying  that  P\  ~  (P0  -  Ps)cxp(-aL),  and  agreement  was  found  to 
within  about  1%  for  zero  chirp. 


3.  Dynamic  Results  for  a  Long  Fiber  Compared  to  Experiment 


Here  we  calculate  Ps  as  a  function  of  Pq  for  a  6-km  single-mode  fiber.  The  FWHM  of  the 
spontaneous  Brillouin  peak  was  measured  to  be  AvB  =  39  MHz  using  an  optical  spectrum 
analyzer  (OSA).  Also,  the  fiber’s  modal  area  was  measured  by  scanning  a  razor  blade  across  the 
beam  in  the  far  field  of  the  exit  face.  The  mode  field  radius  was  then  calculated  to  be  r  = 

4.55  pm  at  the  exit  face. 


The  reflectivity  Ps/Po  is  plotted  in  figure  1  for  both  an  unchirped  and  a  chirped  laser.  It 
approaches  100%  for  large  incident  powers,  illustrating  that  Brillouin  scattering  limits  the 
transmission  through  the  fiber.  The  threshold  (defined  as  the  incident  power  at  which  the 
reflectivity  is  equal  to  1  %)  is  approximately  two  orders  of  magnitude  larger  in  the  chirped  case 
than  in  the  unchirped  case.  For  incident  powers  well  below  threshold,  the  reflectivity  levels  off 
to  a  spontaneous  value  that  is  independent  of  chirp  and  is  in  good  agreement  with  the 
experimentally  measured  value  of  Rq  =  (3.0  ±  0.5)  x  10~6,  with  no  free  parameters. 


incident  laser  power  (W) 


— • —  0  Hz/s 
— le!4  Hz/s 


Figure  1.  Brillouin  reflectivity  versus  incident  power  for  both  an  unchirped  (J3  =  0)  and  a  chirped  (at  />  = 
1014  Hz/s)  laser  pumping  a  long  fiber. 
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To  directly  compare  these  results  to  the  experimental  data,  the  total  backscattered  power  is  next 
plotted  in  figure  2  as  a  function  of  Pq.  This  backscattered  power  is  the  sum  of  the  Brillouin 
Stokes  power  Ps  and  the  Rayleigh  power  PR  =  2.3  x  10~4  Pq.  The  Rayleigh  backscattered  power 
was  determined  from  low-power  measurements  with  an  optical  spectrum  analyzer.  There  is 
reasonable  agreement  between  theory  and  experiment. 


incident  laser  power  (W) 


•  exptfor  0  Hz/s 
- sim  for  0  Hz/s 

•  exptfor  lel4  Hz/s 
- sim  for  le!4  Hz/s 


Figure  2.  Comparison  of  the  simulated  and  experimental  (7)  total  backscattered  power  for  the  unchirped  and 
chirped  laser  sources. 


The  peak  value  of  the  Brillouin  gain  coefficient  is 


So  = 


2  ny1 

'y 

ncp0v\AvBM 


(8) 


—12 

which  equals  6.4  x  10  “m/W  for  the  parameters  given  above.  This  value  is  a  factor  of  2  smaller 
than  what  is  typically  measured  for  fibers  (8),  possibly  because  of  inhomogeneous  broadening  of 
the  spontaneous  Brillouin  peak.  According  to  equation  8,  it  is  only  the  product  g0AvB  that  should 
be  the  same  for  all  polarization-scrambling  silica  fibers  pumped  at  1.55  pm  and  not  their 
individual  two  values. 


The  power  spectral  density  (in  dBm/Hz)  of  the  reflected  Stokes  electric  field  Es  (Op)  was 
computed  as  a  function  of  frequency  for  each  simulation  run,  and  plotted  up  to  the  Nyquist 
limiting  frequency.  Lineshape  functions  were  then  fit  to  these  spectra.  The  resulting  FWHM  are 
plotted  as  the  blue  dots  in  figure  3  for  the  long  fiber  at  an  incident  laser  power  Po  of  1  mW,  well 
below  threshold,  for  four  different  chirps.  At  the  two  lower  chirps  (3  of  1010  and  1012  Hz/s,  the 
lineshapes  were  found  to  be  Lorentzian,  whereas  at  the  two  higher  chirps  of  5  x  10  and 
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1  T 

10  Hz/s  the  lineshapes  were  Gaussian.  The  continuous  red  curve  is  a  plot  of  the  broadening 
expected  in  a  simple  model,  namely 

A  v  =  A  vB  +  nL/3  /  c  (9) 

which  fits  the  simulation  results  to  within  the  error  bars  on  the  fitted  widths. 


log  of  chirp  (Hz/s) 


Figure  3.  FWHM  Avof  the  Brillouin  peak  as  a  function  of  the  chirp  ft  on  a  semilog  scale. 


4.  Dynamic  Results  for  a  Short  Fiber  Compared  to  Other  Models 


Consider  a  shorter  fiber  having  a  length  of  L  =  17.5  m,  core  radius  of  r  =  13.75  pm,  and 
spontaneous  Brillouin  FWHM  of  AvB  =  20  MHz  (corresponding  to  a  peak  Brillouin  gain  of  go  = 
1.2  x  10-11  m/W),  as  might  be  used  to  passively  deliver  a  high-power  laser  beam.  All  other 
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parameters  are  taken  to  be  the  same  as  those  listed  in  sections  2  and  3.  The  reflectivity  is  plotted 
in  figure  4  for  seven  different  values  of  the  chirp  / 3  in  Hz/s. 


- threshold 

lelOHz/s 
1  el 3  Hz/s 
-*-lel4Hz/s 
-»-5el4Hz/s 
-k-  1  el 5  Hz/s 
-*-3el5  Hz/s 
1  e  1 6  Hz/s 


Figure  4.  Brillouin  reflectivity  for  a  laser  beam  having  a  linear  chirp  ft  ranging  between  10 10  and  10*6  Hz/s  that  is 


incident  on  a  short  fiber. 


Well  below  threshold,  the  reflectivity  has  a  constant  value  of  R0  ~  5  x  10’ 10  independent  of  chirp 
and  pump  power.  Jenkins’  heuristic  steady-state  model  (i)  for  the  buildup  of  the  spontaneous 
Stokes  wave  from  thermal  noise  predicts  that  it  should  scale  linearly  with  the  effective 
absorption  length  of  the  fiber,  La  =  [1  —  exp {-aL)\  /  a ,  and  be  independent  of  A  vfe, 


2n2kTy2La 
Mn2  PqU2  A^A 


(10) 


This  equation  implies  Rq  ~  3  x  10  9 ,  which  is  only  a  factor  of  6  larger  than  our  simulation  value. 
Similar  agreement  is  found  for  the  long  fiber  in  section  3 — equation  10  predicts  Rq  ~  9  x  10  6, 
whereas  figure  1  has  a  low-power  reflectivity  of  Rq  ~  2  x  10’ 6. 


The  curves  in  figure  4  cross  the  1%  reflectivity  line  at  the  threshold  incident  laser  powers,  Pth- 
Those  crossing  points  are  plotted  in  figure  5  as  the  blue  dots.  At  the  maximum  chirp  of 
1016  Hz/s,  the  threshold  has  increased  by  a  factor  of  50  compared  to  an  unchirped  pump  source. 
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Figure  5.  Threshold  power  Pttl  as  a  function  of  the  chirp  p  on  a  semilog  scale. 


In  figure  5,  the  threshold  power  scales  linearly  with  the  chirp  ft  above  1014  Hz/s.  One  expects  the 
threshold  to  increase  as  t(/tc,  where  the  fiber  transit  time  is  q  =  nL/c  and  the  time  required  for  the 
pump  laser  to  chirp  out  of  the  Stokes  bandwidth  is  tc  =  Avq/2/3.  Combining  a  resonant 
integration  of  this  adiabatic  increase  (7)  with  the  familiar  factor  of  21  for  the  unchirped  threshold 
value  of  g0PLL/A,  one  obtains 


21A  tfltc 
g()L  tan_1(rf  !tc) 


(ID 


Equation  13  is  plotted  as  the  red  line  in  figure  5.  The  agreement  with  the  time-dependent 
simulations  is  excellent,  with  no  freely  adjustable  parameters.  One  can  easily  verify  by  graphing 
both  sides  that 


tf!tc 

tan-1(ff  !tc) 


(12) 


for  values  of  tf/tc  ranging  from  0  to  6,  thereby  showing  that  equation  9  is  consistent  with  the 
adiabatic  model  on  which  equation  1 1  is  based. 
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5.  Conclusion 


Time-dependent  numerical  simulations  show  that  SBS  in  optical  fibers  can  be  suppressed  by 
linearly  chirping  the  pump  laser  center  frequency.  One  significant  advantage  of  this  method 
compared  to  competing  techniques  that  broaden  the  pump  laser  with  white  noise  is  that  a  narrow 
linewidth  is  maintained  at  any  given  instant,  enabling  coherent  combination  of  the  outputs  of 
several  such  fiber  lasers  to  scale  up  to  higher  powers.  Linear  chirping  enables  a  mismatch  in  the 
length  of  one  such  amplifier  relative  to  another  to  be  compensated  by  shifting  the  frequency  at 
the  input  to  that  particular  fiber  using  a  modulator. 

These  dynamic  simulations  are  in  good  agreement  with  experimental  results  for  a  6-km  fiber  in 
the  laboratory.  The  simulations  should,  therefore,  also  be  valid  for  a  short  17.5-m  fiber  having 
characteristics  appropriate  to  high-power  laser  delivery  systems.  The  time-dependent  computer 
results  for  such  a  short  fiber  are  in  good  agreement  with  a  simpler  adiabatic  model,  thereby 
validating  use  of  that  model  to  predict  the  performance  of  chirped  fibers  without  requiring  the 
demanding  computer  overhead  of  the  full  time-dependent  calculations. 
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